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Abstract 



We study a coarsening model describing the dynamics of interfaces in the one- 
dimensional Allen-Cahn equation. Given a partition of the real line into intervals of 
<*' , length greater than one, the model consists in constantly eliminating the shortest 

^ I interval of the partition by merging it with its two neighbors. We show that the 

CS| ' mean-field equation for the time-dependent distribution of interval lengths can be 

. explicitly solved using a global linearization transformation. This allows us to derive 

I rigorous results on the long-time asymptotics of the solutions. If the average length 

O ' of the intervals is finite, we prove that all distributions approach a uniquely deter- 

r~| ■ mined self-similar solution. We also obtain global stability results for the family of 

"j^ I self-similar profiles which correspond to distributions with infinite expectation. 

. 1 Introduction 

X 

5t ! Consider a domain D C M" which is divided into a large number of sub domains (or 
cells) of different sizes, separated by domain walls, and assume that the system evolves 
in such a way that the larger subdomains grow with time while the smaller ones shrink 
and eventually disappear. In particular, the average size of the cells increases, so that the 
subdivision of D becomes rougher and rougher. Such a coarsening dynamics is observed 
in many physical situations, especially near a phase transition when a system is quenched 
from a homogeneous state into a state of coexisting phases. Typical examples are the 
formation of microstructure in alloy solidification ||LiS61|| and the phase separation in 
lattice spin systems [pe97| , |KBN97|] . Closely related to coarsening is the coagulation (or 



aggregation) process which describes the dynamics of growing and coalescing droplets 
DGY91| , PeR92| , |Vo85|| . In this case, the system consists of a large number of particles of 



different masses which interact by forming clusters. Again, the total mass is preserved, 
so that the average mass per cluster increases with time. 
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Given a coarsening or a coagulation model, the main task is to predict the long-time 
evolution of the size distribution of the cells, or the mass distribution of the clusters. In 
many cases, experiments and numerical calculations show that this behavior is asymp- 
totically self-similar: the system can be described by a single length scale C{t), and the 
distribution approaches the scaling form C{t)~^^{x / £{t)) as t ^ oo. The profile $ and 
the asymptotics of C{t) can sometimes be determined exactly [[NaK86| , BDG94 |. However, 



even in simple situations, it is very difficult to prove that the size distribution actually 
converges to a self-similar profile. 

In this work, we consider a simple coarsening model related to the one-dimensional 
Allen-Cahn equation dtu = d^u + |(u — u^), where a; G M. The equilibria of this system 
are the homogeneous steady states u = ±1, together with the kinks u{x) = ±tanh(x/2) 
which represent domain walls separating regions of different "phases" . If u is any bounded 
solution of this equation, then for t > sufficiently large the graph of u(t,-) will typically 
look like a (countable) family of kinks separated by large intervals on which m ^ ±1. If 
we denote by Xj(t) the position of the j^^ kink and if we assume that Xj+i(t) — Xj(t) ^ 1 
for all J e Z, a rigorous asymptotic analysis shows that xj ~ F{xj+i — xj) —F{xj —Xj-i), 
where F{y) = 24e~^ [|CaP89 |. In other words, the positions of the domain walls behave 



like a system of point particles with short range attractive pair interactions. Thus, on 
an appropriate time scale, only the closest pairs of kinks will really move; in such pairs, 
kinks will attract each other until they eventually annihilate. 

This kink dynamics suggests the following coarsening model |[NaK86|, pGY91|, |GaP92|, 



BDG94 |RuB94| , |BrD95| , |CaPOO|| . Consider a partition of the real line M into a countable 



union of disjoint intervals Ij, with > 1 for all j G Z. In the previous picture, the 

intervals Ij correspond to regions where u is close to ±1. A dynamics on this configuration 
space is defined by iterating the following coarsening step: choose the "smallest" interval 
in the partition, and merge it with its two nearest neighbors. This model clearly mimics 
the dynamics of the domain walls in the one-dimensional Allen-Cahn equation. However, 
proving that the formal procedure described above actually defines a well-posed evolution 
(e.g. for almost all initial configurations) and investigating its statistical properties after 
many coarsening iterations is a non-trivial task, which has not been accomplished so far. 
Instead, the coarsening model has been studied in the mean field approximation, which 
consists in merging the minimal interval not with its true neighbors, but with two intervals 
chosen at random in the configuration {Jjjjg^. This approximation is valid provided 
the lengths of consecutive intervals stay uncorrelated during the coarsening process, an 
assumption that seems reasonable pDG94]| . 



Under this hypothesis, it is possible to write a closed evolution equation for the distri- 
bution f{t,x) (per unit length) of intervals of length a; > 1 at time t [paP92|| . Denoting 
by N(t) = f(t,x) dx the total number of intervals per unit length, and by C{t) the 
length of the smallest interval, the equation reads 

dtf{t, x) = ^(^){(^' -^(^)) ( r ""^'^ f^t, y)f{t, x-y-C{t)) dy - 2/(t, x)N{t) 1 , (1.1) 



N{ty 



'0 



for X > C(t), whereas f(t,x) = for x < C(t) by the definition of C(t). By construc- 
tion, N{t) decreases with time, while the total length of the intervals Jg xf{t,x) dx is 
conserved. 
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We prefer to work with the distribution density p{t,x) = f(t,x)/N(t), which satisfies 
pit, x) = for X < C{t) and the normahzation p(t, x)dx = 1 for all t. The evolution 
equation for p reads 

rx-Cit) 

dtp{t,x) = C{t)p{t,Cit)) p{t,y)p{t,x~y~C{t))dy ioix>C{t). (1.2) 

Jo 

Of course, systems (|1.1| ) and (|1.2|) are equivalent. In particular, once the density pit, x) is 
known, the total number N{t) can be recovered by solving the ordinary differential equa- 
tion N{t) = -2t{t)p{t, C{t))N{t), and the distribution f{t, x) is then given by N{t)p{t, x). 

It is important to note that equations ( |1.1| ), ( |1.2| ) are invariant under reparametriza- 
tions of time. As a consequence, the minimal length is not determined by the initial 
data, but can be prescribed to be an arbitrary (increasing) function of time. In ||CaP92|| , 



the authors define an "intrinsic time" by imposing the relation f(t, C{t))C{t) = 1, which 
means that the number of merging events per unit time is constant. We find it more 
convenient to use the "coarsening time" defined by the simple relation C{t) = t. In other 
words, we choose to parameterize the coarsening process by the length of the smallest 
remaining interval, forgetting about how much physical time elapses between or during 
the merging events. With our choice, equation ( |1.2|) becomes 

PX — t 

dtp{t,x) = p{t,t) p{t,y)p{t,x-y-t)dy ioi x > t. (1.3) 
Jo 

Since we do not allow for intervals of length smaller than 1, we impose our initial condition 
at time t = 1: p(l,x) = Pi(x). 

The aim of this paper is to show that the dynamics of ( |1.3| ) can be completely under- 
stood using a global linearization transformation. As a consequence, we are able to prove 
that solutions of (|1.3| ) satisfying xpit,x) dx < oo approach a non-trivial self-similar 
profile as t — > oo. To achieve this goal, we first rewrite ( |1.3| ) in similarity coordinates by 
setting 

p{t,x) = -r]{\ogt,x/t), or rfir^y) = pie" ^e'y), 

where r = logt > and y = x/t E [1, oo). Then the rescaled density ?7(r, ■) lies in the 
time-independent space 



P 



|r/GLi((l,oo),M+) I ^ r/(y)dy = l}, (1.4) 

which is a closed convex subset of L^((l,oo)). Moreover, (|1.3| ) is transformed into the 
autonomous evolution equation 

drVi'T.y) = dy{yr]{T,y)) +r]{T,l) j 7]{t, z)r]{T,y-z-l)dz for ?/ > 1. (1.5) 



In Section y we show that, for all initial data rjQ G P, ( p..5|) has a unique global solution 
Tj e C0([0, oo), P) with ri{0) = rj^. 

We now define a nonlinear map A/" : P ^ L\^^i[l, oo), R^) by 
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where T is the Fourier transform and 0(z) 
in P, a direct calculation reveals that w(t^ 
drw{T,y) = dy{yw{T,y)). As a consequence, 

w{T,y) = {SrWo)iy) = 



ilogY^f- Vi^r) is a solution of (|T]5|) 
= Nirjij, •)) satisfies the linear equation 



ewQ^ey) if y > 1, 
if2/<l, 



where Wq = Af{rio). It follows that any solution r] G C°([0, cxd), P) of ( [L5|) satisfies 
J\f[rj{T)) = SrJ^{i]o) for all r > 0. In other words, the nonlinear evolution defined by (|1.5|) 
is conjugated (via the map Af) to the linear semigroup (Sr). Thus, the difficulty of solving 
( |1.5D is carried over to the study of the mapping A/" and of its inverse A/""^ = o0~^ ojF. 
Although the properties of these maps are not fully understood, it possible to obtain some 
information on them using the analyticity properties of the Fourier-Laplace transform. 

In Section ^ we investigate the steady states of (^5), which form a one-parameter 



family {rj^jg^^. Here t]l = M ^{^w*), where w*{y) 



,-1 



y ^fell- 



More explicitly, we have 



Vm = {^Vim = tanh - Ei(iO for ^ G 



e 



;i.6) 



where Ei is the exponential integral ||AS72|| . We prove that r^g G P if and only if G (0, 1]. 



Moreover, ?7*(?/) decays exponentially as y ^ oo, while r]g{y) ~ y ^^^^^ if < < 1. In 
particular, r^* is the only steady state for which the average length f^yTil{y)dy is finite. 

Finally, Section |^ is devoted to the convergence results. If the initial data rjo E ¥ 
satisfy y'^rjo G L^((1,cxd)) for some 7 > 3/2 (so that yrio{y) dy < 00), we prove that 
the corresponding solution of ( p..5|) converges exponentially to the steady state r]^: 



?7lj||L2((i,oo)) 



0(e 



-(7-3/2)rN 



for r 



00. 



In terms of the original variables, this shows that the density p{t, x) asymptotically 
approaches the self-similar solution t~^ril{x/t) of (|1.3| ). Moreover, the remainder is 
(9(^^-{7-3/2)^^ so that the convergence is very fast if 7 is large, i.e., the initial data decay 
rapidly at infinity. Similarly, if < < 1 and if 7]q eF satisfies y'^{i]o~^Ve) ^ -^^((l? ^)) 
for some 7 > 6 + 1/2 and some > 0, we prove that the solution of (|1.5| ) with initial data 
rjo converges to the steady state rjg. 

To conclude this section, we briefiy comment on previous results and possible general- 
izations. The mean field equations ( |1.1|) and especially the self-similar solutions ( |1.6|) can 
be found in many physics papers |[NaK8(j| , |DaY9l| , [BPGH |RuB94| , |BrD95|| . The first 
mathematical work is ||CaP92|] , where the authors prove the existence of global solutions 
to ( |1 . 1| ) . They also show that the profile rjl is a positive function (a crucial property 
that is tacitly assumed in the physics literature!) and study its asymptotic behavior as 
y ^ 00. Our main contribution is the introduction of the linearization transformation A/" 
which allows to prove the convergence results. We also extend the analysis of | UaP92 | to 
the equilibria rjg with < 6 < 1. 

The "two-sided" coarsening model discussed in this introduction is clearly not the 
most general system to which our analysis applies. For instance, we can consider the 
"one-sided" variant in which the minimal interval is only merged with one of its neighbors 
JaPOOII. More generally, we can assume that, for j = 1, . . . , A^, the minimal interval has 
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a probability pj of being merged with j of its neighbors, where pi + ■ ■ ■ + = ^- In the 
mean field approximation, this leads to an evolution equation similar to ( |1.3|) , where the 
quadratic convolution in the right-hand side is replaced by a more general convolution 
polynomial. Except for a modified definition of the mapping Af, this extension does not 
affect our analysis in any essential way. Therefore, in the rest of this paper, all results 
will be stated and proved in this general situation. 

Acknowledgments. The authors are grateful for financial support through the French- 
German grant: PROCOPE 00307 TK, Attractors for Extended Systems. 

2 The coarsening equation and its solution 

As is explained in the introduction, we shall study a general coarsening model for which 
the number of intervals involved in each merging event is not necessarily fixed. Instead, 
we allow for some randomness by choosing nonnegative real numbers pi, ■ ■ ■ ,Pn satisfying 
Pi + ■ ■ ■ + Pn = 1, where pj is interpreted as the probability for an interval of minimal 
length to merge with j other intervals. We define the polynomial 

N 

which satisfies Q(l) = 1. The original coarsening model related to the Allen-Cahn equa- 
tion corresponds to the particular case where Q{z) = . 

If p e LH^)> we set 

AT 

Q[p] = Ep^/^*'' (2.1) 

where p*^ = p * p * ■ ■ ■ * p [j factors) and * denotes the convolution product in L^(M). 
In particular, we have Q[p](x) dx = Q{J^ p{x) dx). In what follows, we shall mainly 
use the space P of probability densities defined by (|0|)- Any p G P can be extended to 
the whole real line by setting p{x) = for x < 1. This natural extension, still denoted 
by p, will be used in the sequel without further mention. As an example of this abuse of 
notation, if p G P, we have Q[p] G P and support(Q[p]) C [2, oo). 

The problem we are interested in can now be stated as follows. Given pi G P, find a 
density p : [l,oo)^ satisfying p(l,x) = pi{x) for a; > 1, p(t,x) = for 1 < a; < t, 

and 

dtp{t, x) = p{t, t)Q[p{t, •)] (x-t) ioi x>t>l. (2.2) 

If Q{z) = z"^, the evolution equation ( |2.2|) reduces to ( |1.3| ). 

By assumption, the density p{t, x) is nonzero only in the sector {{t,x) G 1 1 < t < x}, 
where it satisfies ( p.2|) . An important role will be played by the values of p on the 
boundaries of this domain, namely the initial density pi and the trace of p on the diagonal 
X = t, which we denote by a: 

a{t) = p{t,t) fort>l. 
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Any sufficiently smooth solution of ( p.2|) satisfies p(t, ■) G P for all t > 1 provided 
Pi G P. Indeed, it is obvious from ( p.2|) that p stays nonnegative. Moreover, if m{t) = 
Jt° ^) ^ direct calculation shows that 

^m{t) = a{t){Q{m{t)) -1) for t > 1. (2.3) 

Therefore, if m(l) = 1, then m(t) = 1 for all t > 1. 

A very remarkable property of equation (2^) is that it can be explicitly solved using 
Fourier (or Laplace) transform. If p G P, we define 



/oo 
e-«^p(x)dx for^G 



Then p G C°(M, C) satisfies p(0) = 1, |p(OI < 1 for all ^ ^ 0, and p(0 ^ as ^ ^ ±oo. 
Moreover, p is a positive definite function (in the sense of Bochner). Since support(p) C 
[l,oo), the Fourier transform p can be continuously extended to the lower complex half 
plane 

L- = G C I Im^ < 0}. 

This extension (still denoted by p) is analytic in the interior of L~ and satisfies the bound 
|p(OI < e^""^ for all ^ G L". 

Remark. The closely related Laplace transform is defined by 

/oo 
e"^^p(a;)da; for Rep > 0, 

so that p{p) = p(— ip). In the sequel, we prefer using Fourier transform instead of Laplace 
because the inversion formula is more natural. 



Applying Fourier transform to (|2.2| ) and using the fact that convolutions are turned 
into multiplications, we find the equation 

9ip(t,0 = «(t)e-'«*(Q(p(t,0)-l) fort>l, (2.4) 

where a{t) = p{t,t). To solve (P]W), we introduce the nonlinear complex transformation 
defined by 

Remark that (p'{z) = YlT=oiQ(^)]^ ' ^^^at has a power series expansion with non- 
negative coefficients whose radius of convergence is equal to 1. In particular, the map 
: [0, 1) [0, oo) is one-to-one and onto. Let ip = (p^^ be the inverse map, which 
satisfies 

i/j'iw) = 1-Q{tlj{w)), ^(0) = 0. (2.6) 

By construction, ip is analytic in a neighborhood of the real positive axis. In the particular 
case where Q{z) = z"^, one finds 

1 I + z 

(j){z) = - log- and iplw) = tanh(w). 
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Applying the nonlinear transformation simplifies equation ( p.4|) a lot. The function 
w(t,^) = (j){p(t,^)), which is defined at least for Im^ < 0, satisfies the differential equation 

dtw{t, = -a(t)e-'«* for t > 1, 

which has the explicit solution 

w{t,i) = w{l,0 - a(s)e~'^Ms for t > 1 and Im^ < 0. (2.7) 

Remark that |p(t,OI < e*^"^ for all t > 1 and all i G L", because p(t, ■) G P and 
support (p(t, ■)) C [t, oo). Since 0(z) = z + 0{\z\'^) as z ^ 0, it follows that \w{t,^)\ = 
I0(p(i, 0)1 ^ as t ^ oo if Im^ < 0. Thus, taking the limit t ^ oo in (^.7)), we find 
u}(l,0 = /i°°Q;(^) e~^^*dt, which in turn implies 

/oo 
a{s)e-'^'ds for t > 1 and Im^ < 0. (2.8) 

This formula has a very nice interpretation. Let Af be the nonlinear transformation 
defined (at least formally) by 

jsf = ^ ^ ^ jz or = jr-i o ^ o J^^. (2.9) 

Setting t = 1 in (p^ ), we obtain 0(pi) = S, that is a = N{pi). In other words, the trace 
a(t) = p(t, t) is obtained from the initial density pi(x) = p(l, x) by applying the nonlinear 
map M . Moreover, if U (t) is the linear operator defined for t > 1 by 

(c/(iWM = i,,>„»w = { (2-10) 

then ( |2.8| ) reads w{t, ■) = (p{p(t, ■)) = jF(f/(t)a), which means J\f{p(t, ■)) = U{t)a. There- 
fore, the solution of ( p.2|) satisfies 

M{p{t, ■)) = U{t)M{pi) for t > 1. (2.11) 

This shows that the dynamics of the nonlinear system ( |2.2| ) is conjugated via the nonlinear 
mapping M to the linear evolution U . Since M{pi) is the trace function defined by 
a{t) = pit, t), it is very natural that the evolution of a is obtained just by cutting off the 
history in [1, t). 

It is not difficult to show that the map A/" is well-defined on the space P, cf. ( |1.4D : 

Proposition 2.1 // p G P, then A/'(p) G -^[^^([1, oo), IR+), and the mapping p i— ^ A/'(p) is 
one-to-one. 



Proof. For p G P we construct w = M{p) as follows. Define w : ^ C by ■u;(0 = 
(/)(p(0), where L~ = L~ \ {0} . We recall that p is continuous on L~, analytic in the 
interior of L~, and that |p(OI < 1 for ^ 7^ 0. Since (j) is analytic in the unit disk of 
C, it follows that w is continuous on L~ and analytic in the interior of L~. Moreover, 
< 0(|p(OI) — 0(6^™^)) hence = C(e^™^) as Im^ — >• — oo. These properties 
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imply (see |[Sch66|] , Ch. VIII) that w is the Fourier transform of a uniquely determined 
distribution w G V'(R) with support in [1, oo). 

The injectivity of A/" follows from the facts that the mapping <^ : {2; | |2;| < 1} ^ C is 
locally injective (as 0'(^) = ^/ i^~Q{^)) 0) that : [0, 1) — *• M is globally injective 
(as (p'{s) > 1 for s G [0, 1)). If A/'(pi) = Af{p2), then, by the above, we have 4>{Pi{0) = 
0(^2(0) ^ ^ This proves pi(— ip) = p2(— ip) for p > 0, as pj(— ip) G [0,1). By 
continuity of and local invertibility we obtain pi = p2 on L~, and hence pi = p2- 

To prove w = Mi^p) G L\^^{\l,oo)), choose any e > and consider the distribution 
We : X ^ e~^^w{x). It belongs to iS'(R) (the space of tempered distributions) and its 
Fourier transform satisfies 

^^.(0 = w{i-ie) = 0(p(e-i£)) for Im^ < 0. 
Now we observe that p(^— le) = ^s{i)) where Pei^x) = e~'^^p{x). Since UpeUli < < 1, 
the series ^"^r^P*^ converges in L^(M) to some function G ^^"'^([1, 00), IR+). (Here 

we use the crucial fact that (f)^''\0) > for all k G N.) By construction, 

= E ^^(P.(0)' = 0(p(e-i^)) = ^^.(0 for Ime < 0, 

k=l 

giving We = We G -^^"'^((1, oo),M_,_), and hence w : x e^^We{x) lies in Ljq^([1, oo),M_,_). ■ 
Remarks. 



1. Under the assumptions of Proposition one has that w = M{p) G 5'(R), i.e., w 



is a tempered distribution. In fact, there exists a constant C > such that |iB(^)| 



l0(p(O)l — C*rnax{l, — log 1^1} for 7^ 0, see the proof of Proposition ^?T] below. This 
means that the singularity of w{^) at ^ = is (not worse than) logarithmic. 
2. More information on M can be extracted from the proof of Proposition |2.1| . For 
instance, if p G P, then M{p){x) = p(x) for a. a. x G (l,n+l), where 

n = mm{j e {1, . . . , N} \ pj > O] > 1 (2.12) 

is the largest integer such that |Q(-2)| = 0(|-2|") as z 0. Indeed, in view of (^.5]), one 
has (f){z) = z + C(|;z|"+^) as z ^ 0. It follows that 

~ 0(fc)(O) 

k=n+l 

where the second term in the right-hand side is supported in the interval [n+1, 00). Thus 
We = Pe a.e. in [l,n-|-l], which proves the claim. Similarly, using the observation that 
support(p*'^) C [k,oo), it is easy to show that, if p : [1, cxo) R_(_ is continuous, so is 

^fip)■ 

The formula ( p.llD is very nice, but does not provide an effective method for solving the 



Cauchy problem associated with (|2.2| ). Indeed, Proposition ^TT] does not give a sufficient 
characterization of the set A/'(P), which is also the domain of J\f~^. It is not even clear a 
priori that this set in left invariant by the linear evolution U(t). For this reason, we shall 
use standard PDE techniques to prove existence of solutions to ( p.2|) in the next section. 
But the representation ( p.ll| ) will be very useful to find self-similar solutions of ( |2.2| ) in 



Section ||, and to study their stability in Section 



8 



3 The Cauchy problem for the rescaled system 



The evolution equation (|2.2|) is not autonomous, and it is defined on the time- dependent 
domain {x G M+ \ x >t}. These drawbacks are eliminated if we rescale the density p{t, x) 
by setting 

p{t,x) = jri{\ogt,x/t) foTx>t>l, (3.1) 

or equivalently 

r]iT, y) = eXe^ e^y) for r > 0, y > 1. (3.2) 

In what follows, we denote by r = logt and y = x/t the new time and space coordinates. 
The rescaled density rji^r, ■) now belongs to the fixed space P defined in ( |1.4| ). Moreover, 
it satisfies the autonomous evolution equation 

drV{r,y) = dy{yr]{T,y)) + l3{r)Q[7]{r,-)]{y-l) for y > 1, (3.3) 

where (3{t) = r]{T, 1) is the new trace which relates to a{t) via (3{t) = e'^a(e'^). The initial 
condition for ( |3.3|) is ri{0,y) = rio{y), where ?7o = Pi ^ 

The nonlinearity in ( p. 3D has the form /5(r)TiQ[r7(r)], where Ti : P — P is the shift 
operator defined by 

iT.m-{"Y> ll^^, (3.4) 

In particular, for all G P, the support of TiQ[r7] is contained in [2, oo), or even in 
[n+l,oo), where n > 1 is defined in ( |2.12|) . Thus, any solution of ( p.3|) satisfies the 
linear equation drf] = dy{yr]) in the strip {(t, y) | r > 0, 1 < y < 2}. It follows that 
ri{T,y) = e^~'^°rj{TQ,e^~'^°y) for all r > tq > and all ?/ > 1 such that e^~'^°y < 2. Setting 
?/ = 1, we obtain the important relation 

(3{t) = e^-^° r7(ro,e^-^«) for < t-tq < log 2, (3.5) 

which means that the trace /3{t) for r G [tq, tq+ log 2] can be determined from the solution 
//(tq, ■). This formula will be useful to define the trace (3 properly when the solution Tjij,-) 
of ( |3.3| ) is not continuous. For instance, if ri{T, ■) G P for all r > and if (3 satisfies 
then/3GLi„,([0,oo),M+). 

The main purpose of this section is to show that ( ^.3] ) defines a well-posed evolution 
in the space P. To do this, we consider the associated integral equation 

v{r) = SrVo+ [ f3{s)Sr-sTiQ[v{s)]ds forr>0, (3.6) 
Jo 

where (5't-)t>o is the linear semigroup on P defined by 

(S,,)(y) = { ^'t' Z<l. 

To formulate our convergence results in Section ^, we shall need some weighted spaces 
which we now introduce. For p G [1, oo) and 7 > 0, we denote by the function space 

LP = {we Lf„,([l,oo),M) I \\w\\,^^ < 00}, (3.8) 



9 



where 

. i/p 

\w\\pn = Ib^^lU" = ( / {y^\w{y)\ydy 

When 7 = 0, we simply write instead of Lq and instead of ||w||p^o- Remark that 
L'P if and only if 7 > 1 — 1/p (when p > 1) or 7 > (when p = 1). In what follows, 

we shall often restrict ourselves to such values of p, 7. 

We first give a few basic estimates on the semigroup (Sr) and the nonlinear ity Q acting 
on LP. 

Lemma 3.1 Let p G [l,cx3) and 7 > 0. Then (|3.?1 ) defines a strongly continuous semi- 
group (S't-)t->o in L^, and 

\\SrV\\p,'y<e-^^''-'-''^''^Mp,y, (3.9) 

for all T] E LP and all t > 0. Moreover, equality holds in (|3.9|) if and only if rj^y) = for 
almost all y G [1, e"^]. 



Lemma 3.2 Let Q be the nonlinear map defined by (|2.1|). 

a) IfriE L\ then Q[ri] G L^ and \\Q[v]\\i < QiWvh)- Ifv,V^ L\ then 

mv]-Q[v]\\i<Q'ir)\\v-v\\i, 

where r = max{||r7||i, ||^||i}. Finally, if rj Ef, then Q[ri] E P. 

b) Let p E [1, 00) and 7 > 1 — If rj E L^, then Q[ri] E L^, and there exists C > 
(independent of rj) such that 

||TiQM|U,<CQ'(||r^||i)||r/|U,. (3.10) 

Ifr],f] E LP and R = max{||?7||p^^, ||?^||p,7}, then 

\\TMv]-TMv]\\prf < CQ'(i?)h-r/||,,,. 

Proof. Estimate ( p.9|) is a straightforward calculation, and the proof of Lemma 3^ will 
be outlined in Appendix 0. ■ 

We are now ready to state the main result of this section: 

Theorem 3.3 For any rjo E L^((l,oo),M) with \\rio\\i < 1, equations (|3l6|), ( |3.5| ) have a 
unique global solution rj E C^{[0,oo),L^), which satisfies \\f]{T)\\i < 1 for all r > 0. In 
addition, 

IJif'noEF, then r]{r) E P for all r > 0; 

^) "^f Vo ^ for some p > 1 and some 7 > 1 — 1/p, then rj E C°([0, 00), LP). 
Proof. Fix r]Q E Bi, where Bi = {ri E L^ \ \\r]\\i < 1}. Setting tq = in p.5| ), we obtain 

/5(r) = e>o(e^) for < r < log2. (3.11) 
The first step is to show that (|3.6|) , (|3.11| ) have a unique solution rj E C°([0, log2], L^). 
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Let q = Q'{1) > 1, and let T = (log2)/m, where m G N* is sufficiently large so that, 
for all A; = 1, . . . , m, 

/ e>o(e^)|ds< -. (3.12) 

J{k-1)T Q 

We introduce the Banach space X = C°([0, T], L^) equipped with the norm 

\\v\\x= sup \\vir)\\i. 

0<T<T 

Let B = {i] E X \ \\i]\\x < 1}, and let F : X \—>- X he the nonlinear map defined by 

{F[r]]){T) = SrVo+ [ mSr-sTiQ[v{s)]ds for < r < T, 
Jo 



where (3{s) is given by ( 3.11 ). We claim that F{B) C B and that F is a strict contraction 
in B. Indeed: 

a) Assume that rj E B. Using Lemmas |3Tl| and |3.2| , we find, for all r G [0, T], 

||(FW)(r)||i < \\SrVo\\i+ riPmSr-sT^Qlvmids 

Jo 

/OC pT 
e^\r]oie^y)\dy + |/3(s)| ||TiQ[r/(s)]||ids 



\Vo{y)\dy+ / e^|r/o(e^)|||Q[r7(.)]||ids (3.13) 

^0 

< / \Vo{y)\dy + QiM\x) [ \vo{y)\dy < 1, 



since Q{\\r]\\x) < <5(1) = 1 and ||?7o||i < 1. This shows that F{B) C B. 
b) If r],fie B, then for all r G [0, T], 



mr)]){r)-{F[f,]){r)\U < / ms)\\\Sr-s{TMv{s)]-TMv{s)])hds 

e^\Vo{en\\mv{s)]-Q[fi{smids 



< / e^\vo{e')\Q'{l)Ms)-v{s)\Uds 
Jo 

< e''|?7o(e^)|dsj||r7-r/||x. 



In view of ( ^.12| ), this shows that F is a strict contraction in B. 

Let G X be the unique fixed point of F in the ball B. Then rj satisfies (|3.6| ), and 
using Gronwall's lemma it is readily verified that rj is in fact the unique solution of (|3.6|) 
in the whole space X = C°([0, T], L^). Repeating the same argument m times (where m 
is such that ( p.l2|) holds), we conclude that equations ( p.6|) , ( p.ll|) have a unique solution 
7] G C°([0, log2], L-*^), which satisfies ||^7(t)||i < 1 for all r G [0,log2]. Moreover, it is clear 
that ( |3.5|) holds for all tq G [0,log2] and almost all r G [ro,log2]. 
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For T G [0,log2], let E^- : Bi ^ Bi be the nonlinear map defined by Sr^o = vi'^)^ 
where ?7(r) is the solution of (|3.6|) we have just constructed. Then it is easy to verify that 
'Eti+t2 = '^Ti ° for < Ti + T2 < log 2. It follows that the family (S,-) can be extended 
to a continuous semiflow {Er)T>o- By construction, if rjo G Bi and if we set r/(r) = ErTjo 
for all r > 0, then rj e C°([0, cxo), L^) is the unique solution of (|3.6|), (|3.5|), and rj^r) G -Bi 



for all r > 0. This proves the first part of Theorem 3.3 



Assume now that r]o G P. Keeping the same notations as above, we define 

B = {rieX\ riir) G P for all r G [0,T]}. 

In particular, i? is a closed subset of B, as P is closed in Bi G . If rj G B, it is clear 
that (F[r7])(r) G oo), M+) for all r G [0,T], and that all inequalities in (|3.13|) can 

be replaced by equalities. Thus F{B) C B, hence the solution rj G C°([0, oo), L^) of (|3.6|) 
satisfies r7(r) G P for all r G [0,T]. Proceeding as above, we then show that ?7(r) G P for 
all r G [0, log 2], hence for all r > 0. This proves assertion 1) in Theorem p.3[ 

Finally, assume that r/o G L?' for somep > 1 and some 7 > 1 — 1/p, and that ||?7o||i < 1- 



Using Lemmas |3.1| , p.2| and a fixed point argument as before, it is straightforward to show 
that the solution r] G C°{[0, 00), L^) of (pD satisfies r] G C\[0,T],Lp) for some T > 
(depending on rjo). Let 

T* = sup{T >0\r]e C\[0, T], L^)} G (0, 00]. 

We claim that T* = 00. Indeed, assume on the contrary that < T* < 00. Since 
lh(r)||i < 1 for all r > 0, it follows from (|3j), ([3J0| ) that 



\\v{r)\\p,,<\\Vo\\p,, + Cq \(3{s)Ms%,,ds for < r < T*. 

^0 

Using Gronwall's lemma and the fact that /? G Ll^^{[0, 00)), we deduce that ||?7(T)||p,-y < C 
for all T G [0,T*). In view of ( p.6|) , ( p.lO|) , this in turn implies that ^^(r) has a limit in 
as r y T*, giving rj G C°([0, T*], L^). Since we have a local existence result in L^, we 
conclude that G C°([0, T], L^) for some T > T*, which contradicts the definition of T*. 



This proves assertion 2) in Theorem |3l3 . 



The nonlinear map A/" introduced in the previous section can also be used to linearize 
]3D. Indeed, the Fourier transforms of p and t] are related via p(t, ^) = ?7(logt,t^), 
so that ( |3.1| ) is just a rescaling of the Fourier variable ^. As is clear from ( |2.9| ), this 
transformation commutes with the action of A/". Thus, if p is a solution of (|2.2| ) with 
initial data pi and if r] is the corresponding solution of ( p.3|) given by ( p.2|) , it follows from 
( CT ) that 

^Ar(r/(logt, Ojlx/t) = A/'(?7o)(a;) for x > t > 1, (3.14) 
where rjo = pi. Setting r = logt and y = x/t, we obtain the representation formula 

Af{7]{T)) = Sr^iVo) forr>0, (3.15) 

where (Sr) is the linear semigroup ( p.7| ). The last result of this section shows that this 
formula is indeed correct: 
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Proposition 3.4 Let rjQ G P, and let f] G C°([0, oo),P) be the solution of ( p.6|) given by 
Theorem gj. Then M{vir)) = 5^7V(r/o) for all r > 0. 



Proof. We establish the formula by returning to the unsealed variables {t, x) and by 
showing that the formal steps of Section |^ can be made rigorous for the solutions of 
(pl). Define p : [1,00)^ ^ M+ by p{t,x) = \r]{\ogt, x/t) if x > t > 1 and p{t,x) = if 
1 < X <t. Then p G C°([l, 00), P), and rescaling ( |3.6| ) we find 

p{t) = U{t) (^pi + ^ a{s)TM[p{s)] ds^ for t > 1, (3.16) 

where pi = rjo E P, a(t) = j/5(logt), [/(t) is the linear operator (|2.10| ), and Tg is the shift 
operator defined as in ( |3.4| ). To simplify the notation, we set f{s,x) = (TsQ[p(s)])(x). 
Then / G (^"([l, 00), L^), so that (s, x) ^— a(s)/(s, x) G Li^qj,([1, cx)), L^). By construction, 
the trace a satisfies the identity 

a{t) = pi{t) + a{s)f{s,t)(\s for set > 1. 

We now apply the Fourier transform to (|3.16|) . For any ^ G L~ and any t > 1, we find 

Pi(x)e-'«Mx + y^ ly a(s)/(s,x)ds|e-'«^'dx. 

Since pi G P, the first term in the right-hand side is absolutely continuous with respect 
to t, and 

/oo 
pi(a;)e-'«^da; = -pi(t)e-^«* for aet > 1. 

The second term can be decomposed as hi(t,^) — h2(t,C,), where 

hi{t,0 = J^^ {Jla{s)f{s,x)ds'je-'^^dx = J\{s)f{s,0ds, 

h2{t,0 = ^ 1^ a(s)/(s,x)ds}e-'«Ma;. 
Clearly, hi(t,C,) is absolutely continuous with respect to t, and 

dMt,0 = = «(t)e-'«*Q(p(t,0) for aet > 1. 

Next, since f{s,x) = for x < s, we have a{s)f{s,x) ds = a{s)f{s,x) ds, and 
this expression is a locally integrable function of x. It follows that h2(t,^) is absolutely 
continuous with respect to t, and 

dth2{t, = e-'«* ^ a{s)f{s, t) ds for get > 1. 

Summarizing, we have shown that, for any ^ G L~, the Fourier transform p(t,^) is abso- 
lutely continuous with respect to t and satisfies 

dtp{t,0 = -e''^'(piit)+ f a{s)f{s,t)ds)+ait)e-^^'Q{pit,0) 



a{t) e-'«* (Q(p(t, 0) - 1) for set > 1. 
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This gives (|2.4|). Now, proceeding exactly as in Section we deduce that ( |2.8| ) holds 
for all if: > 1 if Im^ < 0, and this in turn is equivalent to (|2.11| ). Finally, using the 
transformation ( |3.14| ) we obtain ( p.l5| ). ■ 



4 Properties of the steady states 

This section is devoted to the time-independent solutions of ( p.3|) in the space P defined 

by (0). 

Definition. We say that rjo eF is a. steady state of ( p^) if the solution t] G C°([0, oo), P) 
of ( p.6|) given by Theorem p.3| satisfies t]{t) = rjo for all r > 0. 

The steady states of (|3.3|) will also be called "equilibria" or "stationary solutions". 



Lemma 4.1 IfrjoE^isa steady state of (|3.3|) , there exists f3 >0 such that ?7o(l/) = P/u 
for almost all y G [1, 2]. 

Proof. If rji^r) = r/o, ( |3.6|) implies that ?7o(y) = ^^Voi^^v) for all t ^ [0,log2] and 
cey G [l,2e~'^], because the nonlinearity in ( p. 61 ) vanishes identically for such values of 
T,y. We define F : x t— > rjQ{y)dy > and obtain 

F{x+y) = F{x) + F{y) for x,?/ > and x+y < log 2. 

Since F is continuous, we conclude that F{x) = f3x for some (3 > 0. Differentiating 
imphes P = e^?7o(e^) for sex G [0,log2] which gives the desired result. ■ 

Let rjo E F he a steady state. Since rjo coincides almost everywhere in [1, 2] with a 
continuous function, the constant (3 in Lemma can be identified with rio{l). Clearly, 
the trace function defined by ( |3.5| ) satisfies /3(r) = (3 for all r > 0. In particular, the 
integral equation (^]6|) reduces to 

Vo = SrVo + f3[ S,T,Q[Vo]ds for r > 0. (4.1) 
Jo 

From 7]q eF we now conclude that (3 > 0. 

On the other hand, if r^o G P and w = J\f{r)o), it follows from Propositions |2.1| and 



that r]o is a steady state if and only if SrW = w for all r > 0. In view of {p^, this is the 
case if and only if there exists /?' G M such that w = P'w*, where 



But since w{y) = r^od/) for asy G [1, 2] (see Remark 2 after Proposition pj]) , we necessarily 
have P' = P = r7o(l)- 



Finally, since equilibria are time-independent solutions of (|3.3| ), we certainly expect 
them to solve the ordinary differential equation 

{yvy{y) + P{TMv])iy) = o forz/>i, vii) = P- (4.3) 

Remark that the initial value f3 also appears as a parameter in front of the nonlinear term. 



It is not difficult to show that (|4.3| ) has global solutions: 
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Lemma 4.2 For any /9 G M, equation ( |4.3| ) has a unique global solution 77 : [1, oo 



Proof. For any G N, let Ik = [kn + 1, {k+l)n + 1], where n G is defined in ( p. 12 ) 



For any rj G Ll^^{[l, 00), M), the nonhnear term {TiQ[7]]){y) only depends on the values of 
r]{z) for z < y — n. In particular, {TiQ[ri]){y) = for y < n + 1, so that any solution of 
( [4.3| ) satisfies 77(1/) = [3 /y ioi y G /q = [1, n + 1]. Using this information, one can compute 
(TiQ[?7])(y) explicitly for ?/ G h = [n + l,2n + 1], and then solve ( [4.3|) on this interval 
to determine 77(7/) for y E h. By construction, is smooth on both /q and /i, but rj 
has a discontinuity of order n at ?/ = n + 1, in the sense that the derivatives rj^^^y) are 
continuous for A; = 0, . . . , n — 1, whereas T]^^\y) has different limits to the left and to the 
right at y = n + 1 (if /? 7^ 0). Iterating this procedure, we find that ( [4.3| ) has a unique 
global solution 77 G C""^'^([l, 00), R), which satisfies r] G C°°{Ik) for all A; G N. ■ 



The following result shows that equilibria of ( |3.3| ) indeed correspond to solutions of 
the differential equation 



Proposition 4.3 IfrjQ^F and (3 > 0, the following assertions are equivalent: 

a) r/o is a steady state of (p.3|) with 770(1) = P- 

b) ?7o coincides almost everywhere with the solution of ( |4.3| ) . 

c) A/'(77o) = 

Proof. We already proved that a) ■<=^ c). If 7/0 G P is a steady state with 770(1) = P, it 
follows from ( [4.1| ) that 

^^^^^^ + ^ r^.T.QMd. = 0, 
^ ^ Jo 

for all r > 0. Using (|3.7| ), it is not difficult to verify that the first term converges to (yrjo)' 
in ^'{{1, 00)) as r — >• 0, while the second one tends to /?TiQ[7/o] in L^{{1, 00)). This shows 
that (after modification on a set of measure zero) r/o is absolutely continuous on (1, 00) 
and satisfies the differential equation (|4.3| ) for almost all y > 1. It follows easily that rjo 
is the solution of ( |4.3| ) in the sense of Lemma [1.2| . Thus a) =^ b) . 

Conversely, assume that tiq E F satisfies ( [4 .31) . Applying the semi-group Sr to ([4 .31) 
and integrating over r, we immediately obtain ( [4.1D , which implies that rjo is a steady 
state. This proves that b) ^ a). ■ 

The main goal of this section is to determine for which values of /5 > the solution 7/ 



of (|4.3| ) belongs to P. Our strategy is to use the characterization c) in Proposition |473 . 
Therefore, we are led to study the image of Pw* under the map J\f~^, and this requires 
very precise information on the complex transformations ( |2.5| ) and ( |2.6| ). The following 
quantities, related to the polynomial Q{z), will play an important role in the sequel: 



, = O'(l) > 1 and « = exp (^ - ^-^) d.) < 1. 



(4,4) 



Lemma 4.4 Let 

$(z) = 1 - e-«*(^) for \z\ < 1, 
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where is defined in ( p.5|) . Then $ can be extended analytically to a neighborhood of the 
real positive axis R_|_. This extension satisfies $(-2) > and ^'{z) > for all z > 0. 
Moreover, $(0) = 0, $'(0) = q, $(1) = 1, $'(1) = k, and $(z) ^ R as z ^ 00, where 

R = l + exp( Tf-^ i— -)dz- ^dzV (4.5) 



l-Qiz)J 72 l-Qiz 

Note that R= 00 if Q{z) = z and 1 < R < 00 otherwise. 

Proof. Since the polynomial 1 ~ Q{z) has the unique real positive root z = 1, which is a 
simple root because Q'{1) = g 7^ 0, it is clear that the function 

= exp ^ - j dt J = ^ forN<l, 

can be extended to an analytic map in a neighborhood of the real positive axis M+. 
Moreover, x(0) = 1, = ^! and using z—1 = exp{— shows that {z—l)xiz) — > 

R—1 for z — > 00, where R is defined in (|4.5| ). Since $(2;) = 1 — {l—z)x{z), we conclude 
that the function $ has the desired properties. In particular, 

<l>'(z) = qx{z)- ^ 



l-Q(z)' 

so that ^'{z) > for all z>0. ■ 

It follows from Lemma [4.4| that the map $ : [0, 00) —* [0, R) is one-to-one and onto. 
Let ^ = $-1 : [0, R) [0, 00) be the inverse map. Then ^(0) = 0, ^'(0) = ^(1) = 1, 
^'{l) = and \E''(m) > for all u G [0,R). By construction, 

^{u) =i:(^~\og{l-u)^ forO<M<l. (4.6) 

Lemma 4.5 The function \E' : [0,R) — > [0, 00) is absolutely monotone, i.e. '^^''\u) > 
for all k E N and all u G [0, R). In particular, ^ can be extended to an analytic function 
on the disc \u\ < R, and there exist nonnegative coefficients {'^k)km* such that 

00 

\]>(u) = ^ M/^u'' for \u\ < R. 

k=l 

Proof. Since \& = we already know that \& is analytic in a neighborhood of [0, R). 
We first show by induction that, for all n G N*, there exists a polynomial P„ such that 

^(")(^) = forO<n<L (4.7) 

g" (1— -u)" 



Indeed, differentiating (|4.6| ) and using ( p^ , we obtain 



^^(^) = ^ ^(^y ) forO<M<l. (4.^ 
g(l-M) 
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Thus ( [4.7|) holds for n = 1 with Pi{z) = 1 — Q{z). On the other hand, differentiating 
( ^TTD and using (^ISD , we find, for < m < 1, 

i^n+i(^(u)) 



with P, 



n+l\ 



P:,{z){l-Q{z))+nqP^{z) 



(4.9) 



Therefore, ( [4 .71) is estabhshed. 

We next show that, for all n G N*, there exists a polynomial Rn{z) with nonnegative 
coefEcients such that 

P„(^) = (l-Q(z))(l-^)"-ii?„(^). (4.10) 

Obviously, ( [4.10|) holds for n = 1 with Ri{z) = 1. Combining (|4.9|) and ( |4.10| ), we obtain 
the recursion relation 

Rn+l{z) = A,{Z)R'^{Z) + A2{z)Rr,{z) + {u-l) A-^iz) Rn{z) , 

where the coefficient functions Aj are given by 



Aoiz) 



1 - Q{z) 

I- z 
q - Q\z) 



N 



N 



1 



1 



Z 



1 



z 



i-1 



— z 

N j-1 



(1 



j=2 k=l 



Because of pj > all ^1,^2, A3 are polynomials (in z) with nonnegative coefficients. 
Thus, the same property holds for i?„ by induction over n. 

Since < ^(m) < 1 and < (5(^(m)) < 1 for all u G (0, 1), it follows from (F^) and 
( CT ) that ^(")(m) > for all n e N and all u G (0, 1), hence also for u G [0, 1]. By a 
classical result of Bernstein (see |[Fe71|| , Section VII. 2), the power series 



00 ^ 
Y] where = tt^^^Ho) > 0, 



k=l 



(4.11) 



converges absolutely and uniformly for \u\ < 1, and defines an analytic continuation of 
\1/ to the unit disk. Moreover, if i?i > 1 denotes the radius of convergence of the series 
( [4.11| ), it is well-known (see for instance |[Ru87|| , exercise 16.1) that the analytic function 
defined by (|4.11|) has a singularity at u = Ri. Since ^'(m) — 00 as u y R, it follows that 
R = Ri. This concludes the proof. ■ 

Example. To conclude this study of the mappings $ and we give an explicit example 
of a nonlinearity Q for which these functions can be calculated explicitly. Let Q{z) = 
{l—a)z + az"^, where a G [0, 1]. The value a = 1 corresponds to the coarsening equation 
( |1.3|) , while a = is a particular case of a model studied in [PaPOO |. Then q = 1+a = 
R = 1 + 1/a, and 



1 + a 



log 



1 + az 
l-z'' 



1 - e"^"' 
1 + ae-i"" ' 
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The auxiliary functions $, \1/ are: 



^ , , (1+0)2; u 

1 + az l+a — au 



We are now ready to state and prove the main result of this section. 



Theorem 4.6 (Steady states of (|37^)) 

Fix 9 > and let rjg : [1, 00) ^ M. be the solution of ( [4 .31 ) with (3 = 6 /q. Then 

a) G P if and only if0<9<l. 

b) If 9 E (0, 1], eF is positive and strictly decreasing, so that yrjliy) as 00. 

c) IfQ<9 <l, then 

y^oo K 1 [1—9) 

where T is the Gamma function and 7e = — r'(l) ^ 0.577216 is Euler's constant. 

d) If 9 = 1, then 

yr^l{y)dy = —. (4.13) 



K 

Moreover, if degQ > 1, there exists A > such that 



For Q{z) = z we have 



lim = -A. (4.14) 



lim = -1. (4.15) 

y^oo y log y 



Remark. It follows from Theorem [4.6| and Proposition [4.3| that ( p.3|) has a unique steady 
state ril such that yr]*{y) dy < 00. 



Proof. We first show that rj^ElP if 0<9<1. According to Proposition [4.3| , it 
is sufficient to prove that there exists an element of P (still denoted by rjg) such that 
^f{7];) = {9/q)w\ Since U'^ = J^-^ o ^ o T and ^(w) = ^(l-e""'^) by (H), this 
relation is equivalent to 

r^ = vl;(l-e-^-*), (4.I6) 
where rjl = Trig and w* = J-'w*. In view of 



w*{0= / dy = E,{iO, (4.17) 

'1 2/ 



where Ei is the exponential integral, see ||AS72||. It is well-known that 



Ei{z) = — log z — ■je + x{^) foi' |arg2;| < TT, (4-18) 

where x : C ^ C is an entire function with x(0) = and x'(0) = 1- Thus, w* is analytic 
in the interior of L~, where L~ = G C | Im^ < 0}. Moreover, Re('u;*(^)) 00 as 
^ — >■ within L~. 
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In Appendix 0, we prove that |1 - e-^^*(?)| < 1 for all ^ G \ {0} and all 9 E (0, 1], 
see also Figure |A.1| . From Lemma we also know that \1/ is analytic in the disk of radius 
R> 1 centered at the origin. Therefore, the map rj^ defined by ( [4.16| ) is continuous over 
L~ (with ril{0) = 1) and analytic in the interior of L~. In addition, since < \u\ 

whenever \u\ < 1, we have the bound 

IrTeiOl < |l-e-'"*^^^| < 2e\w*iO\ for ^ G \ {0}. 

In particular, = 0(e^™^) as Im,^ ^ — oo. By the Paley- Wiener Theorem (see for 

instance ||Ru87|| ), we conclude that rjg = J-'^^rjl G -^^^((1, oo)). 

To prove that rfl is nonnegative, we argue as in [|CaP92|| . Consider the Lapl ace 
transform rjl = Cr]g, which satisfies Tjg{p) = rjl{—ip). As is well-known (see ||I'e71|| , 
Section XIII. 4), positivity of rjg is equivalent to complete monotonicity of rjg, namely 
{-l)%'^''\p) > for all A; G N and p > 0. Recall that 

^(p) = ^(I_e-^'2*(P)) = ^(i_e-^Ei(p)^ ^^-^g^ 
We apply Lemma below with 

. . / (0,1) - M, r (0,oo) ^ (0,1), 



By Lemma [4.5| , /i is completely monotone, thus it remains to show that g[ is completely 
monotone. Observe that g'l = f2 ° 92, where /2 : M ^ M is defined by f2{w) = 9e~^ and 
(72:(0,oo)^Mby 

g2{p) = eE,{p) - iog(-E;(p)) = eEi{p)+p + hgp. 

Clearly, /2 is completely monotone, thus (again by Lemma |4.71 ) it remains to prove that 
5^2 is completely monotone. This follows from the representation 

g'^{p) = -9— + 1 + 1 = 1 + (1-^)1 + 9 C e-'^ds. 
P P P Jo 

Thus, we have shown that 77^ G L^((l, cxd), M.,.). Since rj{p) ^ 1 as p \ 0, we conclude 
that ?7g G and 'r]g{y)dy = 1, i.e., T]g G P. 

Now, fix 6* > 1 and assume that rjg G P, where rjg is the solution of ( [4.3|) with f3 = 9/q. 
According to Proposition Mirjl) = {9/q)w*, so that ([4.19|) holds. Thus, in view of 
( [4.18| ), the Laplace transform of rjg satisfies 



Ve 



(p) = ^ (^1 - / e^^^E-'^^P))) = 1 - e^^E + 0{p^+^) for p\0. (4.20) 



Since > 1, it follows that yTjliy) dy = —{rig)'{0) = 0, which clearly contradicts the 
hypothesis r]g G P. This proves a). 

Next, fix 6* G (0, 1]. To prove that rjg is strictly decreasing, it is sufficient to show that 
rig{y) > for all y >l, since y{ril)'{y) +ril{y) < by ( [4.3|) . Assume on the contrary that 
there exists yo > 1 such that rig{yo) = and rig{y) > for 1 < y < yo- It is clear that 
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Ho > n+1, where n is defined in (|2.12|) . Thus, {TiQ[r)]){yo) > 0, hence rjg' (i/q) < by 
( [4.3|) , which contradicts the fact that rig G P. This proves b). 

Assume now that < < 1. In Appendix |B], we prove that the hmit in the left-hand 
side of (|4.12| ) exists. Let L{9) denote this hmit, and let 



He{y) 



Clearly, y'Hgiy) ^ Li9)/9 as y 
' 'He{p)- 



rig{x) dx for y > I. 



cxD. Thus, the Laplace transform of Hg satisfies 



P 



P 



P 



- Hg -]dt^T(l- 



.m 



as p \ 0. 



Since r]*g{p) = e-P-pHg{p) = 1 - p^T{l-e)L{e)/e + o{p^) asp \ 0, it follows from (|420|) 
that r(l-^)L(^)/^ = e^^^/K. This proves ( CT )- 

Finally, let 6' = L Then ([4.19|) , (|4.2CI| ) show that the Laplace transform rjl is analytic 
in the half-plane {pG C| Rep> —A}, where A > is the unique real root of the equation 
1 — e"^^'^"'*'^ = R (ii Q(z) = z, then R = oo, hence also A = oo.) In particular, r]l(y) 
decays exponentially as y — oo, and 



yr]l{y)dy = — . 



K 



If degQ > 1, then A < oo, and the arguments given in | CaP92|] (in the particular case 
Q{z) = z"^) show that ([4.14|) holds. If Q{z) = z, then A = oo and ril{y) = p{y — l)/y, where 
p : [0, oo) is the Dickmann function studied in [|CaPOO|| . From the asymptotics of 

p given there, we deduce that ( [4.15|) holds. This concludes the proof. ■ 



The following lemma was used in the proof of Theorem ^.6| . For its proof see ||I'e71 
Section XIII.4. 

Definition. Let / C M be an open interval, and let / G C°°(/, R). The function / is 
called completely monotone if {—l)''f^''\x) > for all a; G / and all k eN. 

Lemma 4.7 (Composition lemma) 

Let /, J C M be open intervals. If f : J ^ W is completely monotone and g : I J has a 
derivative g' which is completely monotone, then f o g : I ^ is completely monotone. 

Remarks. The (generalized) steady states rjg with 9 > 1 will not be studied in this paper, 
because they do not lie in our function space P. We just mention here a few properties 
that can established using the techniques developed in the proof of Theorem |4.6| . There 
exists a critical value 6^* G (1, oo] such that 

1) If 1 < ^ < 9^., then rjg G L^((l, oo), M) and r]g{y) dy = 1. However, r]g is not a 
positive function. In particular, \\rig\\i > 1, so that rjg does not belong to the unit baU of 

where existence of global solutions is known from Theorem O. 

2) If ^ > 9^, then r/* ^ L\{l,oo),R). 

Moreover, 6*=,, = oo if Q{z) = z, whereas < oo if degQ > 1. In the particular case where 
Q{z) = z'^, one has 9^ ~ 3.24826. These statements are illustrated in Figure |0| . 
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Figure 4.1: The steady state ?/| of the coarsening equation (3.3) with Q{z) = z"^ is represented for 



four values of the parameter 9. The first two graphs {0 = 0.5 and 9 = 1.0) iUustrate the conclusions of 



Theorem 4.6, and the other two {6 = 2.0 and 6 — 3.5) the remarks after Lemma 4.7. The pictures were 
produced using the explicit formula (|1.6|) and a FFT routine to compute the Fourier transforms. 



5 Global convergence results 



In this final section, we use the exphcit representation formula (|3.15|) to study the long- 
time behavior of the solutions of ( p.3| ). In particular, we obtain global stability results for 
the steady states rj^ with < 9 < 1. 

Since the nonlinear map Af, which allows us to linearize (p.3|) , has a simple expression 
in Fourier variables, it is convenient to use L^-based function spaces instead of the L^- 
based function spaces which are more natural for the existence theory. Our basic space 
will be 

P^ = Fr]Ll for 7 > 0, 



where P is defined in (|1.4| ) and in ( pl8|) . Remark that P.y is a closed subspace of if 
7 > 1/2, since ^ L^. 

The image of under the Fourier transform JF can be characterized completely. Let 
Mj be the space of all functions z : L~ — > C satisfying the following three conditions: 

(i) z is analytic in the interior of L~, 

(ii) for each ^2 < 0, the map ^1 t— > z{^i+iC,2) lies in the Sobolev space H'^i 
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-6 



z\\a-y = supe ^^\\z{-+i^2)\\h'<{r) < oo. (5.1) 
6<o 



(In ( ^.1|) , the supremum over ^2 < is always attained at ^2 = 0.) 

Then 77 G if and only if 77 = J^rj G (when 7 = 0, this is just the Paley- Wiener 
theorem, see [[Ru87|| ; the general case follows using the Fourier characterization of the 



Sobolev space H'^(M.).) Moreover, the map r] WvWml = \\v\\h-i{r) is a norm on 
which is equivalent to ||?7||2,7- If in addition r] G P-y, then 77(0) = 1 and ^1 1— > 77(^1) is a 
positive definite function on M (in the sense of Bochner). Furthermore, \ri{^)\ < 1 for all 
e G L- \ {0}. 

Assume now that t] E P-y for some 7 > 1/2, and let w = N'i'i]), namely wl^) = (j){r]{^)). 
Since is analytic in the unit disk, it is clear that w is analytic in the interior of L~. 
Moreover, the fact that (f){z) = z + O^zl"^) as z ^ guarantees that w{C,) has the same 
decay properties as rf{^) as 00. However, since 77(0) = 1 and since (f){z) has a 

singularity at z = 1, we see that w{^) necessarily has a singularity at ^ = 0. This is 
the reason why the nonlinear transformation Af does not map into itself. To handle 
this difficulty, our strategy is to subtract from w{^) a suitable function with the same 
singularity at ,^ = and whose inverse Fourier transform is explicitly known. 

If 7 > 3/2, a natural candidate for this counter-term is ^w*(0 = 0(^(0); where 77* 
is the unique steady state of ( |3.3[ ) that belongs to L^, see Theorem ^^6] . We recall that 
w* is defined in 



Proposition 5.1 Let 7 > 3/2 and t] G Pj. Then X{r]) = ^w* + d with d G L'^^^- 

Proof. We first show that M{ri) G = L^((l,oo)). As explained above, it is sufficient 
to prove that w = (pijj) satisfies (|5.1| ) with 7 = 0. Choose a > sufficiently small so 
that \(t){z)\ < 2\z\ for all 2; G C with \z\ < a. Since rj G L^, there exists 6 < such that 
1^(01 — <^ whenever Im^ <b. Thus 

supe-^2||w(-+i6)||i2 < 2supe-«2||77(-+i^2)||L2 < 00. 

On the other hand, by a variant of the Riemann-Lebesgue lemma, there exists c > such 
that \ri{C,)\ < a for all ^ G L~ with | Re^ > c. Arguing as before, we thus get 



sup / |w(6+i6)| d^i < 00 

b<i2<0 J \^i\>c 



It remains to verify that 



/ l^i^(6+i6)Pd^i < C uniformly in ^2 e [6,0]. (5.2) 
Since w : ^ C is continuous except at the origin, it is sufficient to establish (|5.2|) for 



b, c sufficiently small. Now, as — >■ in L , we have the expansion 
Vi0 = ^-^f^^ + rii0^ with ri(0 = 



if 7>5/2, 
C)(|^|7-i/2) if 3/2<7<5/2, 
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where /i = yr]{y) dy > 1. Using the representation (j){z) = — (1/g) log(l— $(2;)) to- 
gether with the properties of $ hsted in Lemma [4.4| , we thus obtain 

w{0 = --log(l-<l>(r^(0)) = --log(i«:/ie+r2(0), 
q q 

where k = $'(1) and r2{^) satisfies the same bounds as ri(^). This expansion immediately 
imphes ( |5.2| ) if b,c are sufficiently small. Thus, we have shown that A/"(?7) G L^. Since 
obviously ^w* G L^, we deduce that d = M{ri) — ^w* G L^, too. 

To prove that d G L'^_i, it remains to verify that d G H'^^^(R.). Again, by a localization 

argument, it is sufficient to show that d G H^^^{{—c, c)) for some c > sufficiently small. 
If ^ G M, 1^1 < c, we use the representation 

m = mo) - = -^iog( ^;_t^P )- (5.3) 

From (|4.18| ), we know that e''"*'-^-' = i.^e'^'^e"^^'^-*, where x is an entire function vanishing 
at the origin. It follows that d{^) = — (1/g) log(D(^)/^), where D G H'^ {{—c, c)) satisfies 



-D(O) = and -D'(O) = K/ie The claim is now a direct consequence of Lemma ^ 
below. This concludes the proof of Proposition |5.1|. 



Remark. It follows immediately from the proof of Proposition that 

d{0) = c?(?/)d?/ = -(7e - log(K/i)), where /i=/ yri{y)dy. (5.4) 
Ji 1 Ji 



Lemma 5.2 Let 7 > 1, and let I G M. be an open interval containing 0. There exists a 
constant C{1, 7) > such that, for each f G H'^{I) with /(O) = 0, there exists g G H^~^{I) 
such that f{x) = xg{x) for all x E I and 

\\9\\h^-Hi) < C(/,7)||/||h7(/). 

Proof. It is sufficient to prove the claim for / = M (the general case can be reduced to 
this one using a bounded extension operator). If / G H^(M.) and /(O) = 0, the Fourier 
transform / has zero mean and satisfies A^/ G L\R), where A(0 = (1+^^)^/^- Define 
g G L^(R) by its Fourier transform 

f{s)ds = - / f{s)ds for ^GM. 
-00 J( 

Then xg{x) = f{x) for (almost) all x G M. Moreover, since 



r XisV-^\f(s)\ds ife>0, 

A(0^" l?(OI < < I ^ 

I /LA(s)--V(-)|d3 ife<0, 

it follows from Theorem 328 in |PLP59|| that ||A'^^^^||i2 < 2||A'^/||i2, which is the desired 
bound. ■ 



^1- 



We next show that the inverse map M ^ is well-defined in a neighborhood of ^w* in 
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Proposition 5.3 Let 7 > 1/2. There exists e > such that, for all d & L'^ with ||c?||2,7 < 
e, the function N'^^{^w*+d) is well-defined and lies in L^. Moreover, there exists C > 
such that 

\\Af-\lw*+d)-r]lh,,<C\\d\\2,,, 

where r]l = jV^^^^w*). 

Proof. Throughout the proof, we denote by || ■ instead of || ■ W^q the norm on 
defined by (|5.1|) . We first remark that the space Mj is an algebra if 7 > 1/2: there exists 
Ci > such that \\rs\\^ < Ci\\r\\^\\s\\^ for all r, s G H]'. Moreover, as is easy to verify, 
there exists C2 > such that, for all integer k > 1, 

\\r^h < C^2A;^^'''^^||r||^^^(^'^)||r||^('='^\ (5.5) 

where ||r||oo = sup{|r(^)| | C, G L^} < C|kll7 and i{k,'y) = minjn ^N\n>k or n> 7}. 
Assume that d G L'^ for some 7 > 1/2, so that d G Mj. For all ^ G L^, we define 

r(0 = 1 - e-'^*(«) and s(0 = e-^*^^\l-e-'''^^^^). (5.6) 
From ( [4.17| ), ( 4.18|) , it is easy to see that r G H^, and we prove in Appendix ^ that 



< 1. On the other hand, since is an algebra, it is clear that a = 1— e"'^"' G Mj, 
hence s = (1— r)cr G Mj. In addition, if ||(i||2,7 < £ for some e < 1, there exists C3 > 
such that < C3||(i||2,7 < C^e. In particular, ||s||oo < Ce. 

We now fix i?i G (1, -R), where i? > 1 is defined in ([4.5|) , and we assume that e < 1 is 
sufficiently small so that ||s||oo < Ri — 1- We then define r/ G by 

r}=^jj{lw*+d) = ^(l-e-^(i^*+^)) = ^(r+s). 



where ^ is given by ( |4.6| ). From Lemma [4.5| , we know that is analytic in the disk 
{u G C I |m| < R}, with the expansion \l'(u) = J2k>i'^ku''. Since r + s E Mj and 
Ik + -^lloo < Ri < R, it follows from ( [5. 5] ) that the series \E'(r+s) converges in H]', so that 
rf G H]'. By construction, rf = Trj for some rj E L"^ with A/'(?7) = ^w* + d. 

It remains to show that H^] — ^T^IIt — C'||(i||2,7, where rjl = \E'(r), see ( [4.19| ). Using 
\\{r+s)'' - r'^W^ < CiSup{\\k{r+9sf-^\\^ \ 9 G [0,1]} ||s||^ 
and (5.5|), it is straightforward to verify that there exists C4 > such that, for all > 1, 

< C^e+^R'l-^\\s\\^. (5.7) 

Since 



k=l 

it follows that 



I 1 1 



k=l 
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This concludes the proof. ■ 

Remark. Unhke A/", the inverse mapping A/""^ is not positivity preserving. However, if 
in Proposition |5.3| we assume in addition that t] = Af^^{-w*+d) is a positive function, 
then y yrjiy) G L^{{)-, oo)) and 

yr/(?/)d?/ = -e^E-<?rfo^ where rfo = j d{y)dy. (5.8) 
Indeed, on the one hand the Laplace transform rj{p) = r]{—ip) satisfies 

= - f 1 - ^(l-pe^E-^^^^-^^^P))) ^ - e^E-""^" for p \ 0, 
and on the other hand, using rj^y) > 0, we find 

= / yviy) — > / yv{y)(iy = \\v\\i,i forp\o, 

p Jl py Jl 

by Lebesgue's monotone convergence theorem. 



In addition to Lemma 32, the following bounds on the nonlinearity Q[?7] will be used 



to prove our convergence results: 

Lemma 5.4 Fix 7 > 3/2. For any M > 0, there exists C > such that the following 
estimates hold: 

a) For all ri E with ||r7||2,7-i < M and \\ri\\i^i < M , 

\\TiQ[ri\h,,<q\\r]\\2,, + C. (5.9) 

b) If > 2, then for allrj,fj G with ||?7||2,7 < M and ||r/||2,7 < M, 

||TiQM-TiQ[r/]||2,7 < g||r7-77||2,7 + Ch-r/||2,7_i. (5.10) 

Proof. See Appendix p. ■ 

We are now ready to state the main result of this section, which shows that all solutions 
of ( p.3|) in P^ with 7 > 3/2 converge towards the limiting profile rjl. 

Theorem 5.5 Assume that rjo G P^, for some 7 > 3/2, and let rj G (^"([O, 00), P^) be the 
solution of (p.6|) given by Theorem \3. Sj . Then rj is bounded in and there exists C > 
such that 

Mr) - V*ihy~i < Ce~(^~3/2)r 
Moreover, if'~f>2, then 

Mr)-v*ih'y < C(l+r)e-(^-3/2)" /or r > 0. 
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Remarkably, the faster the initial data decay at infinity, the faster the solution con- 
verges to the steady state. For compactly supported data, it should be possible to obtain 
faster decay than exponential. 

Proof. We first prove ( |5.11|) using the representation formula (|3.15|) . By Proposi- 
tion M^rjo) = ^w* + d for some d G L'^_i- Since the semigroup Sr is linear and 
leaves -w* invariant, we have S't-A/'(?7o) = -'w* + Srd. By Lemma ( ^TTD , USt-c/H 2,7-1 < 



e '^/^^'^ 1 1 (i 1 1 2,7-1, so that Srd — > in L"^-^ as r -h> 00. Thus, when r is sufficiently large, 
we can apply Proposition |5.3| which gives 

Mr) - vlh.-i = \\Af-\lw*+Srd)-r]l\\2,,-i < Cie-(^-3/2)r^ 

for some Ci > 0. This estimate holds in fact for all r > with a possibly larger constant 
Ci, which proves ( p. 11 ). Remark that, since rjlr) is nonnegative and ||S'T-(i||i — > 0, it 
follows from ( ^.8|) that ||?7(t)||i,i ^e'^^ as r ^ 00. 

We next show that ||?7(t)||2,7 is uniformly bounded for all r > 0. We already know 
that ||?7(r)||2,7-i and ||^7(t)||i,i remain bounded. Thus, using the integral equation (|3.6|) 
together with the bounds ( |3.9| ) and (|5.9|), we obtain 

h(r)||2,7<e-(^-i/2)n|r/o||2,7+ F Pis)e-^''''/'^^^-^\qUs)h,, + C)ds, 



for some C > 0. Remark that /3(r) = ?7(r, 1) = Af{ri{T))\y=i = ^+e'^d{e'^). Since d G i^^-i, 
it follows that /3(r) = J + ^(r) with e E L\R+). Setting H{t) = e''"'-^/^^^\\r]{T)\\2,-y, we 
thus find 

H{t) <H{0)+ [ (l+g|e(s)|)i/(s)ds + Ce(^-i/2)r fo^ r > 0, 
Jo 

for some C > 0. Since 7 — 1/2 > 1 and e G L^(M+), it follows from Gronwall's lemma 
that H{t) < C2e(''-^/2)r ^^^^ > j^g^^^g h(r)||2,7 < C2 for all r > 0. 

Finally, if 7 > 2, we show that riir) converges to 77* in L^. To do this, we consider the 
integral equation satisfied by r(r) = 77(r) — 77*, namely 

/*7" "1 

r(r) = 5.r(0) + 5._,|£(s)TiQ[r7(s)] + -(TiQ[77*+r(s)] - TiQ[r/*]) } ds. 

In view of ( p.ll| ) and Lemma |5.4| , there exists C3 > such that 

||TiQ[7/(s)]||2,7 < C^3, ||TiQ[r/*+r(s)] -TiQ[r/t]||2,7 < Ib,^ + C3 e-^^-^'/^^^ 

Using Lemma p.l| again, we find that Rir) = ||r(r)||2,7 satisfies the integral inequality 





{Cs\ 


Jo 





Since e(r) = e'^d{e^) with d G L^_i, we have 

1/2, 
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hence there exists C4 > such that 



R{r) < C4e-(^-=^/2)r ^ re-(7-i/2)(-s)^(^)^^_ 

Jo 

Using Gronwall's lemma, we conclude that R{t) < C^^l+r) e^^'^^^^'^^^ for some C5 > 0, 
which is the desired result. ■ 

We now argue that the convergence towards the steady state rjl cannot be faster than 
g-{7-3/2)r ^Yie norm of L^, so that the result of Theorem |5.5| is optimal. To see this, 
we study the linearization of ( p. 3D around rj^. Setting rjlr) = i]l + b{T), we obtain the 
linearized equation drb = Ab, where 

{Ab){y) = (ybYiy) + (iTi(Q'[ry*] * b) + b{mQ[vl]) iv) ■ 

Since we are interested in solutions ri{T) G P-y, we study this operator in the space 

b{y)dy = 

Proposition 5.6 1/7 > 3/2, the operator A on has a = —(7—3/2) in its spectrum. 

Proof. For 6 > 7+ 1/2 we define a Lipschitz function bs G by 

r -1 foTye[l,Ys], 

bs{y)=l -l + {l+{Ys+l)-'){y-Ys) for y e {Ys,Ys + I), 
[ y-^ for y>Ys + l, 

where Y5 is chosen such that bs has mean 0. Note that Ys has a finite limit as 5 \ 7+I/2. 

Our aim is to show that Abs + (7—3/2)65 stays bounded in as 5 \ 7+I/2, while 65 
is unbounded. For this purpose, we compute the asymptotic behavior of Abs as ?/ ^ 00. 
Since rjl decays faster than e~^^ for some A > and since ril{y) dy = 1, we obtain 
{Q'[7]l] * bs)iy) = qbsiv) + 0{y-^-^) as ?/ ^ 00, where q = Q'{1). It follows that 

iAbs){y) = {-6+l)y-'+{y-iy' + 0{y-'-') + 0ie->^y) 
= {-6+2)y-^ + 0{y-^-^) for y ^ 00, 

where the remainder term is uniform in 6 for 6 ~ 7+ 1/2. This implies the estimate 

\\Abs + (7-3/2)65112,7 < {S-i-l/2)\\bs\\2,, + C<2C 

as (5 \ 7+1/2, since \\bs\\2,j ^ l/v/5-7-1/2. This proves the claim. ■ 

To conclude this section, we also give a global stability result for the steady states rjg 
with < 6* < 1. 

Theorem 5.7 Let < 9 < 1 and 9+1/2 < 7 < min{3/2, 26'+l/2}. Assume that the 
initial value 770 G P satisfies rjQ—uril G for some z/ > 0, and let rj G C°([0, cxo), P) be 
the solution of ( p.6|) given by Theorem \3. 5| . Then there exists C > such that 

Mr) - v;\ky-e < Ce-(^-^-i/2)^ for r > 0. (5.12) 
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Remarks. 

1. From ( [4.12| ), we know that rjgiy) ~ y~^^^ as y — > cxd, so that tJq e Ly if and only 
if 7' < 9+1/2. Thus, the assumption 7 > 9+1/2 guarantees that the difference rfQ—urjl 
decays faster than r/g at infinity (otherwise, we could just choose 770 = r]Q, for some 9' < 9, 
in which case rjlr) = rjg, for all r > so that ( p.l2D certainly fails.) For instance, if 
r/o G P n is such that 



where C > and e > 0, then the assumptions of Theorem |5.7| are satisfied for some u 



and 7. On the other hand, the hypothesis 7 < 2^+1/2 ensures that r]g and hence r]Q lie 
in so that //(r) G L'^_g for all r > 0. 

2. Setting formally 9 = 1 in ( [5. 121 ), we recover ( |5.11| ). However, the main difference 



between the two results is the upper bound 7 < 3/2 in Theorem |5.7| which limits the 
decay rate in time of the perturbations. Even for compactly supported perturbations, the 
convergence in ( ^.12|) is not faster than 0{e~^'^), where 6 = min{^, 1—6'}. 



Proof. The proof is quite similar to that of ( ^.11| ), so we just indicate the main differences 
here. Proceeding as in the proof of Proposition p.l|, we first show that A/'(?7o) = -u]*+d 



for some d G L'^_g. In analogy with (|5.3| ), we find 

The crucial point is the behavior of d{^) as ^ ^ 0, which we now analyze. By assumption, 
Vo = ^Ve + C for some ( E L'^^ with ({y) dy = 1—u. From (2]T6|), we have 

= M/(l-e-^-*«)) = 1 - e-^'=^*(«)i/(e-'''^*«)), 

where H : z ^ {l—'^{l—z))/z is analytic in a neighborhood of zero, with i/(0) = 
*'(1) = 1/k. We recall that q-^^^*^^) = (i^fQ^^Q-Sxi^O ^here x is entire, see ( ^^8|) . 
Since 7 < 2^ + 1/2, we deduce that ri : ^ ^ H{e-'^^*^^^) belongs to H^'-%{-c,c)) for 
some c > 0, and that ri(0) = 1/k. Next, we observe that C(0 — " ^^{Oi "where 
r2 G H"'{{—c,c)) and r2(0) = 0. Using the fractional analog of Lemma |5.2| , we conclude 
that r2{C) = e~^'"^^V3(^), where ra G H^^^{{—c,c)) and r3(0) = 0. (Here, we need 
7 < 3/2: if 7 > 3/2, the claim would be false unless r2(0) = 0.) Summarizing, we have 
shown 

VoiO = + m = 1 - e-'''*^^) + r4(0) , 

where r4 G iJ'^^^((— c, c)) and r4(0) = 0. We now apply the inverse map $ = which 
is analytic in a neighborhood of 1 with $'(1) = k. Using the fact that H'^~^ is an algebra, 
we obtain 

^fioiO) = l-e-''''^^\u + r,{0), 

where has the same properties as r4. Since d{C,) = — ^ log(z/+r5(^)), we conclude that 
de H"'-%{-c,c)) with d{0) = -ilogz/. 
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Now, from (|315| ) we have Af{r]{T)) = Sr{^w*+d) = ^w* + Srd, and \\Srd\\2./-e < 
g-(7-6»-i/2)T||^||^^^^^ for all r > 0. Moreover, it is easy to check that Proposition ^.31 and 
its proof remain valid if we replace everywhere w* with 6w*, rjl with rjg, and 7 with 
7' = 7 — ^ < 9+1/2 (as is explained above, this inequality ensures that 77^ G Ly.) Thus, 
we conclude that 



Hr)-V*eh,^-e = \\Ar-\lw*+Srd) ~ r];\\2,^,e < C\\Srd\\2,^^e = 0{e 
as r 00, which is the desired result. 



-(7-e-i/2)7 



A Bounds on the exponential integral 

Let w*{^) = Ei(i^), where Ei(z) = y^^e''''^ dy is the exponential integral. The goal of 
this section is to prove that 



< 1 for ^ G (0, 1] and ^ G L~ \ {0}. 



For 6 = 1 and G M, this property is illustrated in Figure [A.l 
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Figure A.l: (a) The region D C C delimited by the dotted hne contains the curve {€}*{£_) | ^ G K} (sohd 
hne). (b) The curve {1— e~™ | ^ G M} (soUd hne) is contained in the unit disk of C. 



Fix < < 1, and define : L" ^ C by ^^(O) = 1 and Fg{^) = 1 - e'^^*'-^^ for 
^ G L~ \ {0}. Then Fq is continuous on L^, and analytic in the interior of L~. Moreover, 
Fg is uniformly bounded, because |-Fe(OI — 1 + exp(— 6' Re(u;*(^))) and 



/•oo 1 poo 

Reiw*iO)= / -e«^^cos(eiy)dy> / 

Jl y Jn/2 



cos(t) 



dt 



-0.472, 



for all ^ = ^1+1^2 with ^1 G M and ^2 < 0. Finally, since \w*{0\ < ^ 

as ^2 —00, it is clear that |-F6i(OI — > as ^2 —00, uniformly in ^1 G M. Thus, 
by the maximum modulus principle and the Phragmen-Lindelof theory (see e.g. ||Ru87|| , 
Thm. 12.9), it is sufficient to show that (|0|) holds for all ^ G M \ {0}. 
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Let D C Che the open region defined by 

D = {x+iy G C I |?/| < 7r/2 , x + log(2 cos(?/)) > 0}, 

see Figure |A.1| . As is easily verified, w & D implies |1— e"""! < 1. Thus, all we need to 
show is that 6w*{^) G D for all ^ G M \ {0}. Since G -D and D is convex, it is sufficient 
to prove this property for 6 = 1. 

For ^ > we define 

f°° cos(t) , , sin(t) , , 

^(0=/ —^dt and 1/(0=/ — ^dt. (A.2) 



Then w*{0 = Ei(iO = x(0 - wiO for ^ > and w*iO = + for ^ < 0. 

Moreover, \y{C,)\ < vr/2 for all ^ > 0. Thus, it is enough to verify that K{C,) > for all 
^ > 0, where 

^(0 = + log(2 cos(y(0)) for e > 0. 



We first observe that K{C,) > if ^ > is sufficiently small. Indeed, in view of ( [4.18 
we have the expansions 



x{0 = - loge - 7E + 0{e) and yiO = ^ " ^ + Oi^) for ^ \ 0, 

hence K{^) ^ log 2 - 7e > as ^ \ 0. 

We next show that K{C,) > for < ^ < 7r/2. If not, there would exist ^ G 
(0,7r/2] such that K{^) = and K'{^) < 0. In view of (^, K'{^) < if and only if 
sin(0 sin(?/(0) < cos(0 cos(?/(0)- Since < y{^) < n/2, this is equivalent to ^ + y{^) < 
7i/2, or cos(?/(0) > sin(0- Therefore, ^ G (0,7r/2] should satisfy x{^) + log(2sin(0) < 
i^(0 = 0. But this is impossible, because x(0 + log(2sin(0) —>■ log 2 — 7e > as ^ \ 0, 
and 

^ (x(0 + log(2 sin(O)) = ^^^^1^ > for < e < vr/2. 

It remains to show that K{C,) > for ^ > 7r/2. Let 

X = -x{n/2) ^ 0.472 and y = m&x{y{n/2), -yin)} = -y{7v) ^ 0.281. 

(See ||AS72|| for rigorous bounds on x{^),y{^).) Using the definitions ( |A.2| ), it is easy to 
show that |x(OI < x and \y{^)\ < y for all ^ > tt/2. Thus x(0 + log(2 cos(|/(0)) > 
—X + log(2cos(|/)) > for ^ > n/2. This concludes the proof. 



B Asymptotic behavior of the steady states 

Fix 6 G (0, 1), and let v = V*e ■ [1, oo) ^ M be the solution of (gj) with p = O/q. By 
Theorem |4.6| , rj is positive, strictly decreasing, and j'^ Tj{y) dy = 1. The aim of this 
section is to prove that the limit 

L{e) = lim y^^^rjiy) (B.l) 

y^oo 
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exists (and is finite). This is especially easy in the particular case where Q{z) = z. Indeed, 
since q = 1 and Q[ri] =7] in this case, it follows from (|4.3| ) that 

= yri\y) + v{y) + ^^(y-i) > yv\y) + (1+^)^(1/) for y > 2, 

hence y i— >■ y^^^rj^y) is decreasing (and positive) for y >2. In the general situation where 
N = degQ > 1, we need the following estimate: 

Lemma B.l For all y > N , 

.y/N 



>viy)Q'{J^ vix)dxy (B.2) 



Proof. The only property of t] that will be used in this proof is that f] is nonnegative 
and non-increasing. Thus, by linearity and monotonicity, it is sufficient to prove ( |B.2| ) in 
the case where Q{z) = z^ for some J € N, j > 2. For a > 1, we denote 



Dji^) = {{xi, . . . ,Xj) \ 1 < Xi, . . . ,Xj < a} = [l,a 
Sj{a) = {{xi, . . . ,Xj) \ I < Xi < ■ ■ ■ < Xj < a}. 

Then, for y > j, we have 

Q[v]{y) = / r]{xi) ■ . . . ■ ri{xj)6{xi -\ \-Xj-y)d^x 

JD,{y+l-j) 



= j\ / ri{xi) ■ . . . ■ ri{xj)6{xi + ■ ■ ■ + xj — y) d^x, 

where 5 denotes the Dirac measure. To obtain a lower bound, we replace r]{xj) with rjiy) 
in the last integral, and we perform the (trivial) integration over Xj. We obtain 

Wiiy) > U^-)viy) / vi^i) ■ ■■■■ v{xj-i)(^^~^x, 

where Rj-i{y) = {(xi, . . . , Xj_i) | (xi, . . . , y-Xi- . . . -Xj^i) E Sj{y+l-j)}. Now, it 
is straightforward to verify that Rj^i{y) D Sj^i{y/j). Thus 

QW(y) > U^-)v{y) ri{xi) ■ . . . ■ r]{xj-i)d^~^x 

= jviy) / Vi^i) ■ ■■■■ Vi^j-i) 

JDj.^iy/j) 

ay/j X j-1 . rv/j 

r]{x)dxj = viy)Q'yj ?7(a;)dx 

This concludes the proof. 

Combining ( [4.3|) and Lemma we obtain the inequality 



y-i 

6 / \ 

yv'iy) + viy) + -^(y-i) Q'U vi^) dxj < o for y > a^+i. 
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where ri{y—l) may also be replaced by ri{y). It follows that 

y-l 

Since 77 (y) = 1 and g = Q'{1), the right-hand side becomes negative for y sufficiently 
large. This shows that 

supy^+^/'^7]{y) < 00. (B.3) 

Similarly, for y > 

y-l 

^(y'^'viy)) < y'^'vivMy) with e(y) = ^(l - " r/(a;)dx)) > O. (B.4) 

It follows from ( p.3| ) that 1 — f^^'^^^^ vi^) = 0{y^^^^) for y ^ 00, which yields 
e(|/) = and hence 6 G L^((A^+1, 00)). Thus, the differential inequality ( [B^) 

implies that the limit ( [B.l|) exists. 



C Bounds on the nonlinearity 

In this section, we sketch the proofs of Lemmas |3.2| and |5.4|. Without loss of generality, we 



assume here that Q{z) = for some m G (the general case follows by linearity). To 
bound the convolution products, we repeatedly use Young's inequality ||/*5'||p < 
where / G and g & L^. 

Proof of Lemma |3.2| . 

a)lir] e L\ then Q[r/] = r^**" G and \\Q[7]]\\i < M\T = lfv,V^ ^^ then 



- Q[ri] = {rj—r]) * rj * ■ ■ ■ * rj + ... + fj * f] * ■ ■ ■ * (Tj—rj) (C.l) 

(m terms of m factors), hence WQlf]] — Q[^/]||i < mr™^"^||r/ — "/^Hi = Q'('") ||^ — "^Hi- 
b) Assume now that 7] eLP^ L^. For all y > I, 

y^{T^Q[r]]){y) = [ r/(xi) . . . r/(xj(l+xi+ . . . . . . +x^-l/) d'^x. 



where S denotes the Dirac measure. Due to the support property of rj, only the values 
Xi, . . . , Xm > 1 contribute to the integral. For such values, we have the estimate 

. . . +XmV < C{xl+ . . . +xl), (C.2) 

where C > depends on 771,7. Thus, |?/'''(TiQ[?7])| is bounded by a sum of m convolution 
products of the form \y"''i]\ * \r]\*^"^~^\ Taking the norm and using Young's inequality, 
we obtain 

\\TMv]\\p,y < CmMT-'Mp,, = CQ\Mi)\\v\\p,T 

Finally, using the decomposition ( p.l|) and proceeding as above, we find 

\\TMv]-TiQm\p,y < C{mr'^''\\ri-fi\\p^^ + m{m-l)r'^-^R\\ri-fi\\,) 

< C(Q'(r) + i?Q"(r))||r/-r/||,,„ 
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where R = max{||?7||p^^, Hr^Hp^^} and r = max{||r7||i, Since r < R and RQ"{R) < 

CQ'{R), this is the desired result. ■ 

Proof of Lemma |5.4| . The proof follows the same lines, except that ( |C.2| ) is replaced 
with a different estimate, which can be established by induction over m. If 7 > 1 and 
m G N*, there exists C > such that, for all xi, . . . , > 1, 

m 

i=l jjti 

a) If 7 > 3/2 and rj G Py, then 

\\TMv]h-y < m||r/||ri^l|2,7 + Cm||r/||™-1^l|2,7-i 
= Q\l)Mk, + CQ\M\T,i')Mk,-i- 



b) If 7], 7] E P-y, let M = max{||?7||2,-y, 1177112,7}, n = max{||?7||i,i, ||77||i,i} < M, and r = 
maxjllr^lli, Ht^Hi} = 1. Then 

||TiQ[r/] -TiQ[?]]||2,7 < mr"'^^\\ri-i]\\2,y + m{m-l)r"'-^M\\ri-f]\\i 

+ CmrT'-^llr] - ?7||2,7-i + Cm{m-l)r';;'-^M\\r] - fH^i 
= Q'(l)||ry-r/||2,, + Q"(l)M||77-r^||i 

+ CQ'{r,)\\r] - r/||2,7-i + CQ"{n)M\\r] - fjh,i. 

If 7 > 2, the last three terms in the right-hand side can be bounded by CQ'{M) H^]— ^||2,7-i- 
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